Coherent oscillation in a linear quantum system coupled to a thermal bath. 
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We consider the time development of the density matrix for a system coupled to a thermal bath, 
in models that go beyond the standard two-level systems through addition of an energy excitation 
degree of freedom and through the possibility of the replacement of the spin algebra by a more 
complex algebra. We find conditions under which increasing the coupling to the bath above a 
certain level decreases the rate of entropy production, and in which the limiting behavior is a 
dissipationless sinusoidal oscillation that could be interpreted as the synchronization of many modes 
of the uncoupled system. 



Introduction. We explore some extensions of the problem 
of a two level system in a thermal bath |jj and find some 
new qualitative behavior that can be provided by these 
extensions. Two specific examples are the symmetrical 
double well in an environment with temperature on the 
order of the (single-well) level spacing, and the propaga- 
tion of oscillating neutrinos in hot, dense media. In both 
cases we find behavior in which a weak coupling of the 
system to the bath induces some rate of entropy gain, 
but where an increase in the coupling beyond a certain 
value reduces the rate of entropy production, reducing it 
to zero in a limiting case. The limiting behavior can cor- 
respond to the "frozen" conditions of ref. , and many 
other papers, in which the implicit continuous measure- 
ment of a quantity by the system-bath interaction can be 
said to have frozen that quantity in time. However our 
main focus in the present note is on circumstances that 
give instead a kind of dissipationless collective oscillation, 
which gives the appearance of a synchronization of modes 
that had different frequencies in the absence of the cou- 
pling to the bath. We stress, however, that there appears 
to be little direct connection of these results to the large 
literature of synchronized behavior in non-linear systems 
||-|| or of stochastic resonances in forced systems |? | . 

The system itself has two coordinates, which we char- 
acterize as "horizontal" and "vertical". The horizontal 
component is the "spin" coordinate of a two state system, 
or in cases that may be of interest in neutrino physics, a 
"flavor" coordinate for a 3 or 4-state system. We index 
these states with Greek letters. The vertical variable cor- 
responds to a set of energy levels, Bj. The system basis 
states are denoted \E i7 a). The dynamics of the system 
uncoupled to the bath is governed by the Hamiltonian, 



rr-sys _ 



{EtScj + XajiEiME^aHEi,®. (1) 



We take \ a ,p{Ei) <C E^ for all values i, k. We consider 
a density matrix for the system p sys (t) = Trb a th[p s+b (^)], 
where p s+h {t) is the complete density matrix of the 



model. At an initial time t= 
onal in the vertical indices, 



we choose a form diag- 



(* = 0)=X) PaAE i ,t = 0)\E i ,a){E i ,/3\. (2) 



For these initial conditions, one of the results of a master 
equation derivation is that for times t 3> (AE)^ 1 , where 
AE is a measure of the (single- well) energy spacings, the 
operator p sys remains so nearly diagonal in the vertical 
(energy) indices that we can continue to describe the sys- 
tem by a vector in the energy space, p(Ei, t) a ,p, replacing 
(t = 0) by (t) in (||). Henceforth we suppress the indices 
(a,/3) in the horizontal space, in which p{Ei,t) remains 
an operator. For the bath-system coupling, we 
take the dependence on the horizontal (spin-flavor vari- 
ables) to be in a factorized form, H s ^ — (V, where £ 
operates on the horizontal indices and V depends only on 
the vertical coordinates and the coordinates of the bath. 
This coupling is also taken to be weak, in a way that 
allows the problem to be solved, although it can at the 
same time be strong in another sense, for the realization 
of the limits described above. At the end of this note we 
show that for this form of coupling the Bloch equations, 
generalized to include the vertical structure, are of the 
form, 

j t p(E i) t) = -i[X(E i ) ) p(E h t)}+J2Cp( E ht)CnE j) E i ) 

- ^(c 2 p(^,t)+K^,*)c 2 )E r (^^)- ( 3 ) 

All of the elements in (||) are matrices in the horizontal 
space except for the functions T. These functions obey 
the relation T{E j ,E l ) = exp[(Ej - E i )/T]T{E l ,E ] ) in 
consequence of the thermal equilibrium of the bath. For 
the neutrino application, and with the right identifica- 
tions, (||) is the "Quantum Kinetic Equation" derived by 
McKellar and Thomson M. Also Loreti and Balantekin 
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H have used a special case of (||) to discuss the energy 
conserving case of propagation of neutrinos through a 
medium with a (spatially) fluctuating electron density, 
in which case the collision term becomes a pure double 
commutator in the horizontal (flavor) space. For the dou- 
ble well application, we sketch at the end of this note a 
derivation in the context needed for our present results. 

A. Symmetrical double well 

The horizontal space is 2x2, with 03 = 1 for the parti- 
cle to be on the left side and 03 = — 1 for the right side. 
The barrier height is large compared to AE where AE 
is a measure of the energy spacing for the single well. 
The vertical states are the single well energy excitations. 
We choose \ a ,p(Ei) = [a\] a ^g{Ei). The energy split- 
tings g(Ei), which are essentially tunneling rates, have 
strong Ei dependence. For the "spin" dependence of the 
coupling to the bath we take, £ = 1 + 6<t 3 . In the case 
6 = the transition matrix elements of the operator V 
for Ei — > Ej in the left-hand well, considered by itself, 
are the same as those for the right-hand well, as would 
be the case in the dipole approximation of the interac- 
tion with a radiation field. For the case b = 1 only the 
amplitude on the left-hand side of the well interacts with 
the bath. 

The coefficient T(Ei,Ej) is the (single well) rate for a 
state Ei to make a transition to state Ej. When T rs AE, 
where AE is the scale of the vertical splittings, the bath 
induces transitions among the vertical states, and in this 
case our models differ substantively from models that do 
not have the vertical structure. 

B. Neutrinos 

The question is the time dependence of the flavor den- 
sity matrix for neutrinos that oscillate and are also scat- 
tered as they traverse a dense medium, like that of the 
supernova core, or the early universe at a temperature 
greater than 1 MeV. Here the A QiJ g parameters are the 
combination of the mass 2 /(2i?) matrix and the "index of 
refraction" terms, the latter energy independent and of 
order Gw The functions, T{E il Ej) 1 of order G\, are 
the rates of transition for a neutrino to go from state Ei 
to state Ej, as calculated with the matrix £ set equal to 
unity. The choice 6 = would apply to , v T mixing 
in a flavor-blind medium, that is, a medium in which 
and v T have the same interactions with the environment. 
The choice 6=1 corresponds to "active-sterile" scenarios 
in which one of the neutrinos does not interact with the 
environment at all. 

Results for 6 = 0. As an example we use a double square 
well with infinite barriers to the left and right of x = ±8a, 
and with a central barrier of width 2a and height Uq. 
We choose Uq and the particle mass such that twenty 
states are bound with energies less than Uq. The low- 
lying states are thus very nearly the symmetrical and 
anti-symmetrical combinations of separate states in an 
infinite well of width la. The energy splittings, g(Ei) are 
easy to calculate in the limit of small tunneling. For the 
bath we take bosons with u> = vk, equal energy spacing 



2ttv/L where L — > 00, and with creation and annihilation 

operators, a}j,a„. We take H = J2 n a n a n u 'n + Hq VS \ 
and in the system-bath interaction, H s -b = VQ, we take 
the vertical operator V to have off-diagonal matrix ele- 
ments given in dipole expansion as 0, 



= j ^[(al + a n )5j, k + ik n Xj, k (a f n - a n )}. (4) 



where q is a coupling strength. We calculate the off- 
diagonal dipole matrix elements, Xj^ using the infinite 
well wave-functions and define c,- & = qk n Xj f., where the 
later formalism will justify using (near) energy conserva- 
tion to express k n as ±(-Ej — Ek)/v. The coefficients in 
the bath interaction part of M) are then, 



r(Ej,Ek) = 2TTCj^Ck,j 



0{E - E k ) 



9(Ek — Ej) 



,(E j -B h )/T_l 1 



,(Ej-E k )/T 

(5) 



Incorporating these results and solving (^|) for several 
values of q, we show in Fig.l the probability that the 
particle will be on the left side at time t, with the initial 
condition of being on the left and in a thermal distribu- 
tion at the bath temperature, which we take as T = 5 in 
units of the single- well ground state energy. 




FIG. 1. The probability, P, for a particle to be found in 
the left hand well, for a bath temperature T = 5, in units of 
the ground state energy of the single infinite well. The initial 
condition is that the particle is on the left with a thermal 
distribution of energies. The integers (-3,-1,1,3) labeling the 
curves are the Logio of the bath-system coupling strengths in 
arbitrary units. 



* There are diagonal matrix elements of the kx operator as 
well, since the wells are not centered at x = 0. In principle, 
diagonal matrix elements contribute to (^[) in linear order, as 
a modification of the oscillation matrices A, but in the present 
case these contributions cancel when — k is added to k. 
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For very strong coupling, the probability becomes sim- 
ply periodic, as one mode would be in the absence of 
dissipation. This limit can be solved analytically, with 
the result that the oscillation frequency is the thermal 
average of the oscillation frequencies for the separate 
modes |p~0|| . In fig 2 we show the entropy, as defined by 
S = J2j r ^ r [p(Ej)^og(p(Ej))], where the trace is over the 
2x2 space. We note that as the coupling is increased 
from small to moderate values the rate of entropy in- 
crease rises, but that as the coupling becomes very large 
the rate of entropy increase goes to zero. 
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FIG. 2. The entropy S as a function of time for the same 
values of Logio [coupling] used in fig.l. 

We can generalize (||) to apply to a system in which 
the wells contain N fermions with no mutual interactions. 
For example, after making the appropriate modifications 
to the master equation for Fermi statistics, and using the 
same values of parameters as before, but taking TV = 5, 
all initially on the left and in thermal equilibrium, we 
again look at the strong coupling limit, getting results 
very similar to those shown if fig. 1 for the density matrix 
of a single particle. It is surprising that one can prepare 
a state in which, if we looked only at number density, 
a set of Fermi particles could appear to be oscillating, 
all in phase, between the left and the right, as would an 
assemblage of bosons all in an "average energy" state. 

For the case of neutrino mixing in a noisy environ- 
ment, where the vertical space is a continuum, the gener- 
alized Bloch equations, (ph are integro-differential equa- 
tions. But we find that when they are discretized with, 
say, around 30 energy levels, they are easily solved, with 
the results insensitive to the mesh size. The results are 
qualitatively the same as those discussed above. 

There is an amusing generalization of the above con- 
siderations in the extension to a double well with a time 
dependent bias, SHQ ys — e{t)a^. Consider a case in which 
we begin with e > and e ^> for some set of 

low- lying states (but e <C AE). A group of particles 
initially localized on the left, distributed among the low- 
lying states, and not interacting with the medium, will 
stay largely on the left if the bias is maintained at the 
initial value. Now let the bias be slowly lowered until it 



is equal and opposite to the original bias. The particles 
will efficiently be carried to the right hand well. This is 
the analogue of the MSW transition in neutrino physics. 
It is a rather diffuse transition, since the different en- 
ergy levels go through their avoided crossings at different 
times. With a moderate coupling to the bath, the system 
instead settles down rather quickly to half on the right 
and half on the left, with no interesting persistent order. 
However, if we increase the bath coupling sufficiently we 
regain the shift from left to right. The transition is now 
sharper than it was in the case without the bath. De- 
tailed analysis of this type of behavior for the neutrino 
case is given in [[L0| . 

Results for b ^ 0. We briefly summarize some examples. 
In the case plotted in fig. 1, for the strongest coupling, 
we find little change for 6 = 0.001. But b = 0.005 starts 
to poison the oscillation, damping it by a factor of 20% 
in 5 periods. As we move to much larger b, the motion 
becomes more frozen, so that when, e.g., b = 0.5, 98% of 
the probability has remained on the left hand side over a 
time span of five average periods of the 6 = system. A 
popular interpretation of frozen behavior is that the in- 
teractions with the bath are continuously measuring the 
value of (T3 and projecting the wave function onto the 
initial state. Although there is no measurement implied 
by our equations, the results are consistent with this ter- 
minology, since in the case with b = the system-bath 
interaction does not distinguish the right-hand well from 
the left-hand well and it is for this case that we get the 
limiting behavior of undamped simple oscillations with a 
thermal averaged frequency. 

Derivation of the master equation. We focus on the as- 
pects that give rise to the placing of the horizontal ma- 
trices in © and the detailed form of the RHS of §). The 
methods are closely related to those given in |ll|] and |L2j . 
We introduce an interaction picture based on the above 
division of the Hamiltonian. The system-bath coupling 
in this picture is denoted by Hi(t), with vertical matrix 
elements, 



[Hi(t)]j, k = eiHEjV^-iXiE^tyljiEi 



-E k )t 



(6) 



The leading terms in the system+bath density matrix, 
p S j +h , in this interaction picture are generated by the in- 
tegral equation, 



[^ +b (t 2 ),ff / (*2)],H/(*i) • (7) 



p°+ h (t)=p s+h (t = 0) 
t /.ti 

dt\ I dt2 
o Jo 

In the perturbation expansion of (Q) we find a piece 
with an additional power of the quantities Cj k [Et, EX^ 1 ] 
in each higher order. This provides our definition of lead- 
ing terms and the explanation for how a tiny c| k can lead 
to appreciable effects over long periods of time. These 
terms are generated by the terms in which each successive 
new pair of Hi's in the iteration solution puts the state 
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back to the same vertical (and bath) state. These obser- 
vations provide the basis for having omitted, in writing 
(0), all odd terms in Hi, and they also dictate the form 
of the density matrix to be used in solving the equation, 



s+b _ bath 

Pi — P 



Y,Pi{Ei,t)\j)(j\. 



(8) 



1 



exp 



where, 
P bath =[ll( 

i n 

and where, 

Pi{E ,t) = e a (^)V(-Ej,*)e _iA(E3) * 



(£ata n cj n )/Tl (9) 



(10) 



We substitute (||) into (fjj), take the bath trace, and use 
leading order expressions for each term in the double 
commutator, obtaining, for example, 



Tr bath [(j | Hi (h)p s + h (t 2 )Hi (t 2 ) | j) 

k 

x Pl {E k ,ti)e iKEk)tl &~ iX{Ei)t ^ 
Similarly, we have 



S{h - t 2 ) 



(ii) 



Tr 



bath 



01^/(*i)ff/(t2)^ +b (t 2 )|j> = s(h - h) 



x e 



WEiytife-WAh^riE^EME^h). (12) 



The key to separation of the leading terms in the above 
is performing the sum over the modes of the scalar field 
first. After making the transition to the continuum, we 
make replacements of the form, 

-> 2Trf[(E k -E j )]5{t 1 -t 2 )6(E k -E j ). (13) 

where neither the introduced infrared cut-off, w m j n nor 
the residual terms contribute to leading order, as defined 
above. 



Using (11), (|12j) and their counterparts for the two 



other orderings in (6) , doing the t 2 integral (with the 
delta function symmetrically smeared, as can be justified 
by more accurate integrations) and differentiating with 
respect to t gives, 



pi(E 3 ,t) = J2e lXiE3)t <e- lX{Ek)t pi(E k ,t) 



Eh 

^KE k )t r ~i\(E J )t, 



l C,e~ %/ "^ 3,t T{E k ,E 3 ) 
^(e^ E ^C 2 e~ lX{E ^pi{E 3 ,t) 

pi(E J ,t)e lX{E ^ t C 2 e- lX{E ^ t )j2 T ( E ^ E k)- (14) 



Using ( |i"o| ) we now regain (||) . 

To summarize, we have developed the elements of the 
master equation for a bath-coupled system that is a gen- 
eralization in two directions of the system of the ground 
states of a double- well. On the basis of this equation we 
have found new instances in which strong bath-system 
interactions can lead to more, rather than less, orderly 
behavior. We have not begun to address the question 
of identifying a realistic quantum well system in which 
one could make an experimental application of the con- 
clusions. We should emphasize, however, that our limit 
of "large coupling" to the bath is the limit in which 
r 3> A; it is achievable, in principle, for any strength 
of the system-bath coupling (whether the mechanism be 
photons, phonons or collisions with surrounding parti- 
cles), simply by widening the barrier between the wells. 
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